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ABSTRACT 


The  ATOM  '79  (Acoustically  Tracked  Ocean  Mooring)  experiment 
was  designed  to  study  the  natural  background  of  high  frequency, 
high  wavenumber  processes  in  the  upper  ocean  and  to  determine  the 
magnitude  and  effects  of  mooring  motion  on  the  measurements  of 
velocity  and  temperature.  The  mooring  was  deployed  for  about  one 
month  in  the  central  Gulf  of  Mexico  from  December  1979  through 
January  1980. 

It  was  determined  that  translational  mooring  motion  is  not  a 
significant  source  of  error  and  that  torsional  motions  coupled 
with  the  imperfect  current  meter  directional  response  are  a  more 
significant  source  of  error.  Observations  of  the  background  cur¬ 
rents  showed  no  strong  correlation  between  low  frequency  iner¬ 
tial  period  activity  and  high  frequency  activity,  little,  if  any, 
effect  of  the  proximity  of  the  peak  in  the  Brunt-Vaisala  frequen¬ 
cy  profile  on  the  high  frequency  activity,  and  some  significant 
relation  of  wind  forcing  to  sub-mixed  layer  inertial  currents. 
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MEASUREMENTS  OF  HIGH  FREQUENCY,  HIGH  WAVENUMBER  PROCESSES  IN  THE  UPPER  OCEAN: 
THE  ACOUSTICALLY  TRACKED  OCEANOGRAPHIC  MOORING  EXPERIMENT  (ATOM  ’79)  FINAL  REPORT, 

PART  I 


INTRODUCTION 


During  the  past  several  years  an  urgent  need  has  arisen  to  determine  the  clima¬ 
tology  of  the  natural  variability  of  the  upper  ocean  (that  is,  from  0  to  600  m)  on 
time  and  space  scales  that  correspond  to  high  frequency  and  high  wavenumber  internal 
waves  and  turbulence.  Because  the  physical  environment  of  the  upper  ocean  is  harsh, 
past  experimental  programs  were  faced  with  technological  problems  that  severely  lim¬ 
ited  the  accuracy  and  duration  of  the  experiments.  Quoting  Halpern  (1973): 

For  many  years  before  the  1972  SCOR  (Scientific  Committee  on  Oceanic 
Research)  Working  Group  21,  Third  intercomparison  test,  current  measurements 
were  made  beneath  surface  buoys  with  the  realization  that  the  data  contained 
spurious  currents  produced  by  vertical  motions  of  the  instrument  high  frequen¬ 
cy  wave  motions  and  cable  vibrations  (e.g.,  Webster,  1967).  Results  from  the 
1972  SCOR  test  (SCOR  Working  Group  21,  1975;  Gould  and  Sambuco,  1975;  Gould  et 
al.,  1974)  indicated  that  current  measurements  made  at  intermediate  depths 
beneath  surface  buoys  moored  in  deep  water  were  undoubtedly  contaminated  by 
mooring  motion.  Soon  afterwards,  it  was  widely  conjectured  that  erroneous  cur¬ 
rent  measurements  would  be  obtained  at  all  depths  beneath  the  surface  follow¬ 
ing  buoy.  Evidence  contrary  to  this  view  has  since  been  provided  by  Zenk  et 
al.  1978,  Saunders,  P.M.  (1976),  Halpern  et  al.  (1974),  and  Pollard  (1974), 

To  date,  probably  as  a  result  of  this  widely  held  view,  there  have  been  only 
two  major  experiments  directly  assessing  high  frequency,  high  wavenumber  upper  ocean 
phenomena  between  the  surface  and  several  hundred  meters  depth.  These  two  experi¬ 
ments  are  the  MILE  experiment  (Mixed  Layer  Experiment)  and  JASIN  (Joint  Air  Sea 
Interaction  Experiment).  The  lack  of  credible  data  has  inhibited  progress  in  deriv¬ 
ing  a  "climatology"  of  upper  ocean  variability  comparable  to  the  Garrett  and  Munk 
model  for  deep  ocean  internal  waves.  Indeed,  we  should  anticipate  that  future  upper 
ocean  measurements  would  yield  results  that  fundamentally  differ  from  the  deep 
ocean.  From  the  knowledge  of  the  conditions  near  the  ocean  surface,  we  assume  that 
the  high  frequency  and  high  wavenumber  fluctuations  are  not  horizontally  isotropic, 
homogeneous  or  "stationary. ”  Absence  of  these  conditions  is  due  to  the  proximity  of 
the  air-sea  boundary,  vtfiich  is  the  primary  zone  of  thermal  and  kinetic  energy  fluxes 
between  the  ocean  and  the  atmosphere.  Weather  events,  which  are  important  external 
sources  of  short-term  variability,  are  not  generally  distributed  homogeneously,  nor 
do  they  propagate  isotropical  ly  with  space  time  stationarity .  Consequently,  the  sim¬ 
plifying  assumptions  applicable  to  the  deep  ocean  fluctuations  cannot  be  applied  to 
the  upper  ocean. 

Another  complicating  factor  is  the  intensity  of  the  vertical  shear,  which  gen¬ 
erally  has  maxima  in  the  upper  hundreds  of  meters.  The  upper  ocean  mean  vertical 
shear  distribution  is  assumed  to  be  continually  modulated  in  space  and  time  by  sur¬ 
face  weather,  inertial  waves,  internal  tides,  and  the  motions  resulting  from  low 
"mode"  internal  waves.  Theoretical  results  (Landahl  and  Criminali,  1977;  Chimonas, 
1978;  and  Thorpe,  1978)  indicate  that  the  modulation  of  the  high  frequency  and  high 
wavenumber  internal  waves  by  low  frequency  motions  may  lead  to  strong  nonlinear  cou¬ 
pling  of  disparate  scales  of  motion  or  to  the  breakdown  of  the  high  frequency  and 
wavenumber  waves.  Experimental  evidence  of  the  high  frequency  instabilities  and 
sporadic  distributions  has  been  provided  by  Ericksen  (1978)  and  Ko rotayev  and  Pan¬ 
teleyev  (1977a,  1977b).  Ericksen' s  results  were  derived  from  observations  in  the 
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nain  thermoc line,  hut  Korotayev  and  Panteleyev  extracted  their  results  from  experi¬ 
ments  in  the  seasonal  thermoc line.  Both  exoeriments  had  technical  problems  that 
created  a  wide  range  of  valid  but  different  interpretations.  This  was  particularly 
true  of  the  work  of  Korotayev  and  Panteleyev,  which  presented  serious  inconsisten¬ 
cies. 

The  MILE  and  JASIN  experiments  were  directed  toward  c  lari f ication  of  upper 
ocean  dynamics  and  air-sea  exchanges.  Although  both  involved  good  internal  wave  band 
measurements,  neither  was  designed  to  analyze  the  variability  of  fluctuations  at  the 
maximum  of  the  Brunt-Vaisa la  frequency  and  higher.  Indeed,  this  is  borne  out  by  the 
i ntercompari son  tests  between  different  types  of  current  meters  on  different  types 
of  moorings  in  both  of  these  experiments  (Halpern,  Weller,  Briscoe,  Davis,  and 
McCullough,  1981).  However,  some  valuable  insights  about  the  qualitative  features  of 
the  high  frequency  variability  may  be  available.  Resolutions  of  internal  waves  im¬ 
bedded  in  a  finestructure  backgrounu  require  a  high  degree  of  spatial  and  temporal 
resolution  that  was  unattainable  by  standard  instrumentation  and  mooring  techniques. 

The  Acoustically  Tracked  Oceanographic  Mooring  experiment  of  1979  (ATOM  '79)  is 
the  first  phase  of  a  comprehens i ve  design  process  to  obtain  accurate  measurements  of 
the  high  frequency,  high  wavenumber  processes  occurring  in  the  upper  ocean.  By  high 
frequency,  we  mean  processes  which  have  periods  between  one  hour  and  the  minimum 
Brunt-Vai sala  period.  We  define  high  wavenumber  processes  as  those  which  have  wave 
lengths  less  than  about  30  m  and  greater  than  about  1  m  in  the  vertical  and  less 
than  1  km  in  the  horizontal  directions.  The  specific  region  of  interest  (in  depth) 
is  near  the  peak  in  the  Brunt-Vaisala  frequency  profile.  This  typically  occurs  near 
the  base  of  the  mixed  layer  which,  depending  on  geographical  location  and  season, 
may  range  from  the  surface  of  the  ocean  down  to  100  m  or  more.  This  region  is 
particularly  difficult  to  study  with  existing  current  meters.  It  is  too  near  the 
wave  zone  for  the  use  of  subsurface  floats  in  general,  and  at  other  times  the  depth 
of  the  mixed  layer  places  it  in  the  zone  where  measurements  from  surface-moored 
floats  are  inappropriate.  For  the  purposes  of  the  particular  experiment  described  in 
this  report  and  to  minimize  the  problems  involved  with  the  use  of  a  surface  mooring, 
we  chose  to  employ  a  subsurface  mooring  in  conjunction  with  a  deep,  winter  mixed 
layer . 

DESIGN  PHILOSOPHY 

NOPDA  has  been  tasked  with  developing  methods  for  measuring  high  frequency, 
high  wavenumber  variability  in  the  upper  ocean  that  can  be  transferred  to  other 
branches  of  the  Navy  for  use  in  routine  surveys  for  obtaining  climatology  of  this 
variability  on  a  geographical  and  seasonal  basis.  The  primary  purpose  of  the  ATOM 
mooring,  therefore,  is  to  provide  input  to  the  design  procedure  that  will  result  in 
the  development  of  such  standard  observational  methods.  In  one  sense  then,  the  ATOM 
project  may  be  thought  of  as  being  an  engineering  test  of  a  specific  mix  of  observa¬ 
tional  instrumentation.  On  the  other  hand,  because  the  processes  in  the  upper  ocean 
are  so  poorly  known,  it  has  a  very  strong  scientific  component  that  is  required  to 
obtain  information  on  the  structure  of  the  motions  in  this  region  of  the  ocean. 

OBJECTIVES 

There  are  five  major  objectives  to  the  present  study;  three  are  primarily  sci¬ 
entific  in  nature  and  two  are  primarily  technical.  These  objectives  are: 
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Scientific 


1.  determine  the  degree  to  which  high  frequency  variability  In  the  upper 
range  of  the  Internal  gravity  wave  frequency  spectrum  Is  strongly  dependent  on  the 
proximity  of  the  peak  In  the  buoyancy  profile  and  on  the  effect  of  modulation  by  low 
frequency  processes; 

2.  determine  the  degree  to  which  low  frequency  processes  are  governed  by 
local  surface  meteorology; 

3.  measure  the  vertical  coherence  of  the  high  frequency  current  and  tem¬ 
perature  fields; 

Technical 

1.  Investigate  the  contamination  of  the  current  and  temperature  data  in¬ 
duced  by  mooring  motion; 

2.  determine  vrfiether  It  Is  possible  to  remove  a  significant  amount  of  the 
mooring  motion  contamination. 

DATA  REQUIREMENTS 

To  meet  the  objectives  current  and  temperature  data  were  required  at  a  rapid 
sampling  rate,  about  one  sample  per  minute.  The  spatial  sampling  had  to  be  on  the 
order  of  about  5  to  10  m  In  the  vertical  and  cover  a  distance  of  around  100  m  near 
the  maximum  of  the  Brunt-Val sala  frequency  profile.  To  determine  the  local  buoyancy 
profile  and  to  provide  Information  on  the  vertical  wavenumber  structure  of  the  cur¬ 
rent  and  salinity  fields  as  functions  of  depth,  profiles  of  temperature  and  salinity 
were  required  during  the  deployment  and  recovery  of  the  current  meter  mooring.  Si¬ 
multaneous  current  profiles  made  with  the  use  of  expendable  current  profilers  (XCP) 
were  also  made  at  the  time  of  deployment  to  obtain  information  on  the  vertical  wave- 
number  structure  of  the  current  shear  field.  Expendable  bathythermograph  and  air- 
deployed  expendable  bathythermograph  surveys  were  required  to  assess  the  mesoscale 
features  of  the  oceanic  circulation  in  the  vicinity  of  the  mooring  at  the  time  of 
deployment.  Meteorological  data  and  wave  data  from  nearby  weather  data  buoys  were 
also  required  for  determining  the  degree  to  vtfilch  the  local  meteorology  affects  the 
low  frequency  and,  ultimately,  the  high  frequency  processes  occurring  In  the  upper 
layers  of  the  ocean. 

Determining  the  effects  of  mooring  motion  on  the  current  and  temperature  data 
was  a  major  factor  In  the  design  of  the  mooring;  this  requirement  made  it  necessary 
to  track  the  positions  of  the  mooring  elements  as  they  changed  with  time.  Strictly 
speaking.  It  Is  Impossible  to  precisely  characterize  the  time-varying  state  of  a 
single  point  mooring,  vrtilch  by  definition  has  an  Infinite  number  of  degrees  of  free¬ 
dom.  The  best  that  can  be  done  with  a  limited  number  of  Instruments  is  to  record  the 
time-position  history  of  various  parts  of  the  mooring.  For  this  experiment,  we  chose 
to  measure  the  position  of  two  points  on  the  mooring  (one  point  being  above  the 
dense  array  of  current  meters  and  one  below  that  array)  using  acoustic  pulse  and 
Doppler  navigation  techniques;  these  techniques  are  complementary  to  each  other.  The 
depth  of  the  mooring  at  four  positions  on  the  cable  was  determined  by  four  tempera¬ 
ture/pressure  (T/P)  recorders.  The  direction  of  the  mooring  and  tension  at  various 
points  on  the  mooring  was  to  have  been  measured  by  three  force-vector  recorders 
(FVRs)  supplied  by  the  C.S.  Draper  Lab  at  M.I.T.  Torsional  effects  of  the  mooring 
were  estimated  from  the  case  orientations  of  the  various  current  meters.  The  case 
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orientation  of  the  vector  averaging  current  meters  may  be  obtained  every  averaging 
cycle  (these  sample  periods  were  15/1 6ths  and  1 5/8th s  of  a  minute).  The  case  orien¬ 
tation  of  the  acoustic  current  meters  were  sampled  once  every  8  minutes. 

Moorings  such  as  these  are  expensive  due  to  the  bulk  of  subsequent  data  proc¬ 
essing  and  the  cost  of  the  instruments.  It  is  therefore  desirable  to  exercise  dy¬ 
namical  models  of  the  mooring  to  simulate  motion  under  realistic  current  conditions 
prior  to  implantation.  One  objective  of  this  study  was  to  obtain  current  data  at  a 
number  of  points  along  the  mooring  which  could  be  used  as  input  numerical  mooring 
motion  models.  The  validity  of  these  models  could  then  be  tested  against  the  ob¬ 
served  parameters  of  the  ATOM  mooring  configuration.  At  present,  only  one  dynamic 
mooring  motion  model  is  known  to  have  been  qualitatively  compared  with  actual  moor¬ 
ing  notion  data.  This  model  was  developed  at  the  Draper  Lab,  by  Chhabra  (1977).  This 
model  did  not  appear  to  accurately  simulate  mooring  motions  at  the  very  high  fre¬ 
quencies  that  were  important  in  this  experiment. 

The  site  for  the  ATOM  '79  mooring  experiment  was  in  the  center  of  the  Gulf  of 
Mexico  at  approximately  26°N  90°W.  This  location  was  considered  favorable  due  to  the 
historical  absence  of  large  currents,  the  presence  of  a  meteorological  data  buoy, 
and  the  presence  of  a  relatively  flat  bottom,  which  should  not  be  expected  to  act  as 
a  generator  for  internal  waves.  It  was  also  located  far  from  the  continenta  relf 
break,  ana  thus  was  not  expected  to  be  near  a  source  of  internal  tidal  gene  .ion. 
Nevertheless,  this  mooring  was  high  risk. 

The  risk  to  the  mooring  was  predicated  on  the  large  number  of  instrumer  used 
in  the  mooring  and  or,  the  physical  risks  to  that  instrumentation.  A  further  r  i- 
volved  the  use  of  relatively  new  current  meter  instrumentation  that  had  be  ily 
marginally  tested.  The  majority  of  the  current  meters  used  in  this  mooring  were  Neil 
Brown  acoustic  current  meters.  These  were  chosen  on  the  basis  of  linearity  in  speed, 
small  digital  sampling  noise,  and  because  of  relatively  good  horizontal  response 
(obtained  through  the  manufacturer  specifications).  Although  prototypes  of  this  type 
current  meter  had  been  tested  rather  extensively  by  McCullough  and  others 
(McCullough,  1978;  Appel  1  ,  1978),  the  lot  used  by  NORDA  consisted  of  the  first  pro¬ 
duction  batch.  Due  to  the  short  preparation  time  for  the  mooring  experiment,  it  was 
not  possible  to  extensively  subject  these  current  meters  to  severe  environmental 
tests  and  burn  them  in  electronically.  Thus,  it  was  not  known  what  the  data  return 
or  data  quality  would  be.  In  order  to  circumvent  the  possibility  of  a  very  poor  data 
return,  a  number  of  vector  averaging  current  meters  were  used  in  the  vicinity  of  the 
densely  instrumented  part  of  the  mooring. 

The  environment  in  which  the  current  meters  were  to  be  implanted  also  presented 
a  number  of  serious  risk  considerations.  The  first  consideration  was  the  maximum  ex¬ 
pected  current.  The  central  area  of  the  Gulf  of  Mexico,  while  having  a  generally  low 
background  current,  is  known  to  occasionally  be  subject  to  very  strong  current  and 
shear  fields  due  to  either  Loop  Current  extensions  or  eddies  shed  from  the  Loop  Cur¬ 
rent.  It  is  now  known  that  currents  as  high  as  several  knots  may  be  present. 

The  second  consideration  involved  ship  time.  The  mooring  was  designed  to  be 
implanted  in  mid-winter.  This  was  favorable  from  the  point  of  view  of  having  a  deep 
mixed  layer,  but  the  incidence  of  relatively  intense  cold-weather  fronts  over  the 
central  Gulf  of  Mexico  left  a  very  small  time  frame  in  *hich  to  implant  the  meters. 
In  fact,  there  was  a  period  of  only  two  days  wherein  the  mooring  could  have  safely 
been  implanted.  Fortunately  the  clear  calm  weather  persisted  long  enough  for 
deployment.  The  third  major  environmental  risk  factor  to  be  considered  was  human 
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interference.  (The  NOAA  Data  Buoy  Office  has  lost  several  data  buoys  due  to 
fishermen  "accidental ly"  cutting  the  mooring  lines.) 

About  six  months  prior  to  deployment  o.  the  mooring,  a  workshop  was  held  at 
NORDA  to  assess  the  risk  factors  involved.  A  number  of  academic  and  institutional 
contractors  were  invited  to  comment  on  our  mooring  design.  Comments  were  made  and 
responses  were  presented.  In  general,  the  information  given  NORDA  Code  331  by  work¬ 
shop  participants  turned  out  to  be  of  questionable  value,  since  the  experiment  in¬ 
strumentation  and  scientific  objectives  did  not  match  well  with  the  participants' 
experiences,  t^iich  had  been  based  on  somewhat  different  requirements. 

DATA  PROCESSING 

SHIPBOARD  DATA 

Three  types  of  shipboard  survey  data  were  taken:  expendable  bathythermograph 
(X3T),  expendable  current  profiler  data  (XCP),  and  conductivity-temperature-depth 
(CTO).  A  semi-automat ic  XBT  logging  data  system  developed  by  NORDA  Code  350  (R. 
Holland,  1980)  and  modified  by  Mark  Bergin  of  NORDA  Code  331  recorded  the  expendable 
bathythermograph  data.  The  only  data  processing  included  in  this  system  was  the 
transformation  of  the  fall  time  of  the  XBT  into  depth  and  the  application  of  the 
appropriate  resistance  to  temperature  formula  to  obtain  the  temperature  as  a  func¬ 
tion  of  depth.  The  data  from  the  Neil  Brown  CTDs  were  recorded  on  audio  tape  for 
later  shore-based  processing.  At  the  same  time  the  raw  CTD  data  were  plotted  and 
monitored  with  a  Hewlett-Packard  9825A-based  system.  The  expendable  shear  probe  data 
were  simply  recorded  at  sea  for  later  shore-based  processing. 

CTD  DATA  PROCESSING 

The  CTD  data  were  transcribed  from  audio  to  a  nine-track  800  bpi  computer  com¬ 
patible  tape.  The  data  were  then  subsequently  processed  on  the  NAV0CEAN0  UNIVAC  1108 
computer  system.  The  basic  processing  consisted  of  translating  the  data  into  engi¬ 
neering  units,  editing  and  filtering  the  data  to  reduce  the  temperature  and  conduc¬ 
tivity  time  response  mismatch,  and  final  incorporation  of  the  temperature,  conduc¬ 
tivity,  salinity,  and  pressure  data  to  1  decibar  pressure  intervals. 

XBT  DATA  PROCESSING 

The  shore-based  processing  of  the  XBT  data  consisted  primarily  of  producing 
waterfall  plots,  first  differences,  and  high-pass  filtered  profiles.  No  further 
processing  has  been  accomplished  to  date.  The  plots  of  the  XBT  data  are  to  be  found 
in  Saunders,  Green,  and  Bergin  (1980). 

EXPENDABLE  CURRENT  PROFILER  (XCP)  DATA  PROCESSING 

The  XCP  data  were  taken  by  the  Johns  Hopkins  Applied  Physics  Laboratory.  The 
data  were  then  transmitted  to  Dr.  Tom  Sanford  of  the  University  of  Washington,  Ap¬ 
plied  Physics  Laboratory,  for  processing  in  the  manner  described  in  Sanford,  et  al. 
(1978).  The  processed  data  were  transmitted  to  NORDA  on  a  1600  bpi  nine-track  EBCDIC 
card  image  format.  These  records  consisted  of  the  vertical  coordinate  (in  decibars) 
calculated  from  the  drop  time,  the  temperature  (in  degrees  Celsius),  the  east  and 
north  velocity  components  plus  five  auxiliary  variables  (for  engineering  purposes). 
Forty-five  XCP  drops  were  made.  Of  these,  25  were  made  concurrently  with  a  CTO  pro¬ 
file.  The  availability  of  temperature  records  in  both  the  CTD  and  with  the  XCP  pro¬ 
vided  a  means  for  correcting  the  fall  times  of  the  XCPs  using  the  temperature  as  a 
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fiducial  variable.  The  details  of  the  analysis  for  correcting  the  fall  rates  of  tne 
expendable  shear  profilers  are  summarized  in  Green  and  Saunders  (1931)  which  appears 
as  Appendix  A  in  this  report. 

RICHARDSON  NUMBER  CALCULATIONS 

The  25  simultaneous  CTD  and  XCR  casts  provided  both  density  and  velocity  infor¬ 
mation  required  for  the  calculation  of  Richardson  number  profiles.  The  expendable 
shear  profiler  data  were  corrected  for  fall  rate  and  interpolated  to  1  dbar  pressure 
intervals.  The  CTD  and  XCP  records  were  then  low-pass  filtered  with  a  cut-off  of 
about  S  m.  The  Richardson  number  profiles  were  then  computed  from  the  low-passed  1  -i 
XCP  and  CTD  records.  Listings  of  all  the  programs  used  in  the  shipboard  and  onshore 
reduction  of  the  shipboard  data  wi  1 1  be  presented  in  future  reports. 

MOORED  INSTRUMENT  DATA 

Acoustic  Tracking  Data 

The  initial  processing  of  the  acoustic  tracking  data,  both  the  pulse  and 
Doppler  data  for  the  upper  and  lower  positions  of  the  receiving  equipment,  was  oy 
Dr.  Robert  Spindel  of  Woods  Hole  Oceanographic  Institution.  The  lower  unit  at  2 JO  n 
nominal  depth  ran  from  the  time  of  deployment  until  the  transponders  failed  on  Jan¬ 
uary  8,  1980,  at  about  1845Z.  The  upper  unit  also  ran;  but  as  it  did  not  stop  at  the 
end  of  tape  mark,  it  overwrote  the  first  part  of  the  data.  Thus,  the  data  on  the  up¬ 
per  jnit  started  at  1730Z  on  December  28,  1979,  and  ended  at  January  8,  1980,  at 
1845Z.  The  sampling  rate  for  the  pulse  data  was  once  every  half  hour.  The  lower  unit 
'nterrogated  the  transponders  at  15  and  45  minutes  after  the  hour  while  the  upoer 
unit  interrogated  the  transponders  on  the  hour  and  half-hour.  The  initial  position 
data  were  computed  by  Spindel  and  transmitted  to  NORDA  on  a  standard  ASCII  tape. 

During  the  initial  processing  it  was  found  that  vAiere  considerable  numbers  of 
gaps  where  the  data  were  missing  entirely.  These  gaps  were  not  noticed  initially  and 
led  to  the  appearance  of  a  number  of  very  large  spurious  speeds,  which  made  the 
spectrum  appear  to  he  much  flatter  than  it  actually  was.  We  also  learned  that  the 
algorithm  used  to  process  the  data  included  an  assumption  that  the  depth  of  tne 
acoustic  receiver  transmitters  was  constant  at  100  and  200  m,  respectively,  for  the 
upper  and  lower  receiver  transmitter  units.  We  subsequently  obtained  the  raw  travel 
time  information  between  each  of  the  transponders  and  receiver  transmitters  on  the 
mooring  line.  The  depth  of  the  mooring  was  estimated  using  the  temperature/pressure 
recorder  data.  These  data  were  then  used  with  the  travel  time  data  to  calculate  the 
horizontal  and  vertical  positions  of  the  acoustic  receiver  transmitters  on  the 
mooring  line.  The  gaps  in  the  time  series  were  eliminated  through  linear 
interpolation  between  pairs  of  "good"  points. 

The  Doppler  data  were  also  processed  by  Spindel  at  Woods  Hole.  They  were  re¬ 
ceived  and  used  without  any  further  processing  to  attempt  to  correct  errors  in  the 
data . 

TEMPERATURE  PRESSURE  RECORDER  DATA 

At  the  termination  of  the  ATOM  recovery  cruise  the  temperature/pressure  record¬ 
ers  were  shipped  to  the  Charles  Stark  Draper  Laboratory  in  Cambridge,  Massachusetts. 
The  data  were  then  translated  into  engineering  units,  recorded  on  magnetic  tape  and 
mailed  to  N0R0A.  Initial  plots  of  the  raw  data  were  also  produced  at  the  Draper  Lab¬ 
oratory  and  sent  under  separate  cover  to  NORDA.  The  depth  and  temperature  limits  of 
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the  temperature/pressure  recorders  were  preset  prior  to  deployment  so  that  maximum 
resolution  could  be  obtained.  There  were  several  periods  during  the  mooring  when  the 
temperature  pressure  recorder  dropped  below  the  upper  pressure  limit  and  below  the 
lower  temperature  limit  on  the  recorder;  thus,  there  were  no  direct  pressure  meas¬ 
urements  at  the  maximum  depths  to  vhich  the  mooring  was  depressed. 

The  initial  temperature  versus  time  a~H  pressure  versus  time  plots  appeared  at 
first  look  to  be  very  similar.  A  plot  of  temperature  versus  pressure  confirmed  this 
conjecture  and  a  calculation  of  the  correlation  coefficient  yielded  a  value  of 
-0.97. 


This  result  implied  that  it  should  be  possible  to  estimate  the  pressure/time 
history  of  the  mooring  from  the  temperature  record  in  the  gaps  where  the  pressure 
sensor  was  out  of  range.  A  linear  least  squares  fit  of  temperature  to  pressure  was 
used  for  the  extrapolation  formula. 

CURRENT-TEMPERATURE  DATA 

The  basic  data  processing  for  ACMs  and  VACMs  is  the  same.  The  data  were  trans¬ 
cribed  from  the  cassette  tapes  of  each  current  meter  on  to  a  nine-track  (ACM)  or 
seven-track  (VACM)  tape  for  further  processing  on  NAVOCEANO's  UNIVAC  1108  computer. 
Next,  the  data  were  translated  into  engineering  units  and  placed  in  standard  FEB 
(Fast  and  Easy  Binary)  files  (Hal  lock,  1980).  The  current  and  temperature  data  were 
screened  and  edited,  corrections  were  made  in  situations  where  there  were  obviously 
erroneous  data  points  attributable  to  electronic  system  noise.  The  corrections  were 
applied  by  linear  interpolation  between  pairs  of  "good"  data  points.  After  the  data 
were  edited;  the  basic  plots  of  the  current  and  temperature  data  were  produced.  Tem¬ 
perature  and  current  autospectra  were  computed  using  the  standard  NAVDCEANO  spectral 
processing  program.  (Cross-spectra  were  not  computed  at  this  stage  due  to  inappro¬ 
priate  frequency  band  averaging  in  the  NAV0CEAN0  program.)  The  basic  current  statis¬ 
tics  were  also  computed  and  are  plotted  in  Saunders,  Breen,  and  Bergin  (1980).  The 
cross-spectral  analyses  presented  in  this  report  were  computed  using  a  quasi¬ 
ensemble  averaging  technique.  Each  data  record  was  broken  up  into  non-overlapping 
segments  512  points  long.  A  Hanning  window  (Nuttall,  1981)  was  applied  in  the  time 
domain  to  each  segment  and  were  Fourier  transformed.  The  raw  spectral  estimates  were 
computed  by  multiplying  one  transform  by  the  complex  conjugate  of  the  other;  these 
were  then  averaged  across  the  quasi-ensemble  of  estimates  to  produce  a  smoothed  es¬ 
timate  of  the  cross-specta  and  the  auto-spectra  for  each  signal  pair.  The  coherence 
and  phase  were  then  computed  in  the  usual  manner  from  the  auto  and  cross  spectra. 
Other  than  the  application  of  the  Hanning  window  to  reduce  long  distance  frequency 
leakage  by  f‘°  (Nuttall,  1981),  no  frequency  band  averaging  was  applied. 

Estimates  of  the  vertical  velocity  of  the  water  were  made  from  the  temperature 
observations  between  pairs  of  ACMs.  The  vertical  velocity  was  estimated  by  the  equa¬ 
tion  (the  deviation  of  which  is  given  in  Appendix  B) 

w  =.  .  31  /  II  .  _1  d£ 
at  /  az  pg  dt 


This  raw  estimate  was  then  corrected  by  subtracting  the  observed  vertical  velocity 
of  the  mooring  as  obtained  from  the  temperature  pressure  records. 
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An  attempt  was  made  to  estimate  the  vertical  component  of  the  eddy  viscosity. 
This  was  done  by  the  following  method  (derived  in  Appendix  11): 


vr  -  w '  u '  .  vc  -  w_)v_) 

Lu  u,  '  v  vz 


The  calculations  aave  widely  scattered  results  with  vr  ,  varying  from  between 
+103  cmZ/s  to  -10^cm2/s.  u’v 

The  calculations  of  and  '-j-  proved  to  be  inconsistent  with  respect 

to  the  length  of  the  averaging  interval7  and  the  position  of  the  current  meter.  We 

surmise  that  either  the  assumptions  on  whicn  the  estimate  for  the  v£'s  were  de¬ 
rived  were  wrong  or,  more  than  likely,  the  current  meter  spacing  was  too  large  to 

allow  for  good  estimates  of  u2  and  U'  W‘  (the  error  is  probably  in  the  IJ '  compo¬ 

nent).  This  second  possibility  is  deduced  from  the  low  vertical  coherences  of  u  and 
v  discussed  later. 

RESULTS 

MOORING  MOTION  CHARACTERIZATION  EXPERIMENT 

One  of  the  major  goals  of  the  ATOM  experiment  was  to  obtain  information  charac¬ 
terizing  the  motion  of  the  mooring.  The  primary  purpose  of  this  effort  was  to  deter¬ 
mine  the  extent  to  which  the  mooring  motion  introduces  "error"  in  the  current  (and 
temperature)  measurements.  The  secondary  purpose  was  to  obtain  information  on  the 
currents  that  cause  the  mooring  to  move,  as  well  as  the  motion  of  selected  points  on 
the  mooring  that  could  be  used  to  test  the  validity  of  numerical  mooring  motion 
models.  The  tertiary  purpose  was  to  obtain  information  on  the  stresses  in  the  moor¬ 
ing,  as  well  as  the  high  frequency  accelerations  that  could  be  used  for  improved  en¬ 
gineering  design.  The  experiment  succeeded  in  the  first  two  purposes  and  was  unsuc¬ 
cessful  in  the  third  due  to  failure  of  the  force-vector  recorders. 

Mooring  motion  can  affect  current  measurements  by  translational  additive  ef¬ 
fects  and  interaction  with  the  instrument  response  character! sties . 

ACOUSTIC  TRACKING/PRESSURE  MEASUREMENTS 

An  integral  part  of  the  ATOM  design  was  the  capability  to  determine  the  posi¬ 
tion  of  two  points  on  the  mooring  using  acoustic  tracking  methods.  A  combination  of 
a  time-travel /pul sed  system  and  a  Doppler  tracking  system  was  used.  The  operation  of 
the  Doppler  system  is  described  in  Porter,  et  al.  (  1973)  and  Spindel,  et  al  .  (1978). 
The  combined  system  is  described  in  Spindel,  et  al .  (1976).  One  pair  of  pulse  and 
Doppler  tracking  recorders  was  located  just  above  the  dense  current  meter  array,  the 
second  was  below  (Fig.  1).  Three  acoustic  transponders  with  CW  beacons  were  implant¬ 
ed  on  the  sea  floor  after  the  mooring  was  deployed.  The  positions  of  these  beacons 
were  then  determined  by  a  ship  survey.  The  survey  positions  are  shown  in  Figure  2 
(R.  C.  Spindel,  personal  communication,  15  May  1980).  Or.  Spindel  was  responsible 
for  preparing  the  instruments  prior  to  deployment,  conducting  the  ship-board  beacon 
survey,  and  translating  the  collected  data  into  a  format,  usable  by  NORDA. 
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FINAL  CONFIGURATION  FOR  ATOM  79 

(Acoustically  Tracked  Oceanographic  Mooring) 
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Figure  1.  Schematic  diagram  of  ATOM  '79  mooring 


9 


(4465,0,3196.8) 


SURVEY  RESULTS  (IN  METERS) 


Figure  2.  Location  of  bottom  acoustic  beacons 
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DATA  QUALITY 


The  acoustic  positioning  system  proved  somewhat  less  than  totally  reliable.  The 
transponders  (at  least  two  out  of  three)  failed  on  8  January  1980  at  1845Z.  Prior  to 
tins  failure,  the  lower  transponder  unit  (S/N  003)  ran  from  deployment,  while  a  me¬ 
chanical  problem  in  the  upper  unit  (S/N  004)  caused  an  overwriting  of  the  beginning 
of  the  tape,  leaving  data  from  the  period  of  28  December  1979,  1730Z,  to  8  January 
1980,  1845Z.  The  pulse  data  were  sampled  at  a  rate  of  two  samples  per  hour,  while 
the  Doppler  data  were  sampled  at  120  samples  per  hour.  At  any  one  time  there  were 
usually  only  two  working  (i.e.,  with  acceptable  data  returns)  transponders  for  the 
pulse  system.  It  was  not  unusual  to  find  instances  where  no  good  pulse  data  were  re¬ 
turned.  Because  we  did  not,  in  general,  have  three  ranges  from  the  transponders  to 
the  recording  units,  a  third  measurement  was  required  to  determine  the  position  of 
the  unit. 

TEMPERATURE /PRESSURE  (DEPTH)  MEASUREMENTS 

This  third  parameter  was  the  depth  of  the  instrument.  Temperature/ pressure 
(T/P)  recorders  of  C.S.  Draper  Laboratory  design  were  located  at  four  points  on  the 
mooring.  A  T/P  recorder  was  located  near  each  of  the  acoustic  recording  units 

(Fig.  1).  As  a  result  it  was  a  simple  matter  to  determine  the  depth  of  the  acoustic 

units  as  a  function  of  time  as  long  as  the  depth  of  the  T/P  recorder  was  less  than 
that  corresponding  to  the  maximum  pressure  setting.  To  estimate  depths  when  pressure 
exceeded  the  dynamic  ranges  of  the  instrument,  we  used  the  tight  correlation  between 
the  temperature  and  pressure.  As  previously  noted,  the  correlation  coefficient  was 
-0.97.  (A  scatter  plot  of  pressure  temperature  is  shown  in  Fig.  3.)  A  linear  least 
square  regression  equation  was  obtained  between  pressure  and  temperature  at  the 
nearest  acoustic  current  meter  over  the  period  when  the  pressure  recorder  was  in 

range.  The  estimates  of  the  pressure  when  the  recorder  over-ranged  were  made  by 

applying  this  equation  to  the  upper  ACM  temperature. 

POSITION  DATA  PROCESSING 

Using  the  pressure  records  and  extrapolated  values,  the  depth  of  the  acoustic 
units  were  computed  in  the  usual  manner  by  integrating 


dp  =  -p  g  dz. 


The  horizontal  positions  of  the  acoustic  units  were  calculated  using  simple  geometry 
(e.g.,  Spindel,  personal  communication,  15  May  1980).  The  significant  gaps  in  the 
position  data  were  filled  by  linear  interpolation  from  the  nearest  "good"  points. 
The  position  data  were  then  edited  to  remove  spikes.  The  position  data  were  time 
differenced  to  obtain  an  estimate  of  the  mooring's  horizontal  velocity.  Figure  4 
presents  the  time  history  of  the  mooring  velocity,  while  the  position  history  is 
shown  in  Figure  5.  From  this  latter  figure,  we  see  that  the  mooring  moved  about  0.8 
km  in  the  north-south  direction  and  about  0.45  km  in  the  east-west  direction.  A 
stereo  plot  of  the  position  of  the  acoustic  receiver/transponder  (R/T)  is  shown  in 
Figure  6.  The  pulse  uata  velocity  statistics  are  summarized  in  Figures  7  and  8. 
Figure  7  is  a  joint  speed/direction  histogram.  The  pertinent  statistics  are: 
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SCATTER  DIAGRAM  OF  TEMPERATURE  VS.  PRESSURE 
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CORRELATION  COEFFICIENT  =-.97 
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Figure  3.  Scatter  diagram  of  temperature  vs.  pressure  on  upper  temperature/ 
pressure  recorder 
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NORTH-SOUTH  DISPLACEMENT  (KILOMETERS) 
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Figure  5.  Observed  horizontal  position  of  the  ATOM  mooring  at  the  upper  acoustic 
tracking  unit 
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CUMULATIVE  FREQUENCY  RELATIVE  FREQUENCY  RELATIVE  FREQUENCY 


DIRECTION  (DEGREES) 


Figure  8.  Summary  histogram  plots  of  mooring  speeds  and  directions 
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1.  Average  speed:  0.38  cm/sec 

2.  Standard  deviation:  0.26  cm/sec 

3.  Maximum  speed:  3.67  cm/sec* 

4.  Vector  record  speed:  0.04  cm/sec 

The  graphs  of  the  speed  and  direction  histograms  and  the  cumulative  distribution 
function  for  speed  are  shown  in  Figure  8. 

Velocity  variance  spectra  were  also  computed  for  both  the  pulse  and  Doppler 
data.  The  combined  spectra  are  plotted  along  with  a  typical  ACM  current  spectrum  in 
Figure  9.  The  Doppler  spectra  appear  to  agree  quite  well  with  the  pulse  spectra  in 
the  overlap  region.  The  spectral  level  of  the  mooring  velocities  is  about  a  factor 
of  20  below  the  current  spectral  level  until  a  frequency  of  about  3-4  cph.  The 
Doppler  spectrum  then  becomes  quite  flat. 

FLATTENING  OF  THE  DOPPLER  SPECTRUM 

In  order  to  understand  the  cause  of  this  flattening,  we  will  consider  the  ef¬ 
fect  of  least  significant  bit  error  in  the  Doppler  system.  According  to  Spindel  et 
al.  (1976),  the  change  in  a  range  position  over  a  single  estimate  is  given  by 


ARi 


Ki  c 
'  4f, 


where  Ki  is  the  count  from  the  quadrant  detector,  c  is  the  average  sound  speed  and 
fi  is  the  beacon  frequency.  For  the  ATOM  experiment 


c  =  1503.6  m/s 
fj  =  13  kHz 


Thus,  the  least  significant  bit  error,  a,  is  given  by 


A 


C 

4fi 


=  2.89  cm 


Spindel  et  al .  (1976)  states:  "the  cycle  counter  makes  an  independent  phase  measure¬ 
ment  every  0.1  second."  Thus,  the  total  displacement  during  a  time  T  is  given  by 


1 0T  1 0T 

6(T)  =  2  5.  +  2  A 

i=l  1  i=l  1 


♦There  was  only  one  such  observation.  The  next  higher  speed  was  about  1.3  cm/sec 
This  was  probably  an  erroneous  spike  that  was  missed  in  the  editing  process. 
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where  6i  is  the  true  displacement  over  each  0.1  sec,  T  is  the  time  in  seconds,  and 
Ai  is  a  random  variable,  uniformly  distributed  over  -a/2  to  A/2.  The  variance  of  A 
is  then 


a 


2 

6 


IOTA2 

12 


The  corresponding  velocity  variance  is  given  by 

2  _  a62  _  10A2 
°V  T2  "  12T 


If  we  consider  that  this  variance  will  appear  in  the  spectrum  as  white  noise,  evenly 
distributed  over  the  frequency  range  [o,  fn]  (fn  =  Nyquist  frequency  =  1/2T),  then 
the  spectral  level  of  the  Least  Significant  Bit  (LSB)  noise  is  given  by 


SLSB 


10A2 

6 


1.66a2 


The  units  are  in  cm2/s2/Hz.  To  transform  to  the  more  usual  units  of  cm2/ 
/s2/cph,  we  must  divide  by  3600  giving: 


SLSB  (f)  =  4.611  x  10-4A2 

=  3.85  x  10a-3  cm2/s2/cph 


This  value  is  plotted  as  a  horizontal  line,  labeled  LSB,  in  Figure  9.  It  is  clear 
that  this  lies  very  close  to  the  flat  portion  of  the  Doppler  spectrum  and  is  proba¬ 
bly  the  cause  of  this  flattening.  (There  is  no  evidence  to  support  or  reject  the  hy¬ 
pothesis  that  the  mooring  motion  is  significantly  below  the  currents  in  the  frequen¬ 
cy  band  3-30  cph.) 

ORIENTATION  EFFECTS 

No  current  meter  has  a  perfectly  linear  response  under  all  conditions  of  flow. 
This  correct  statement,  but  unfortunate  fact,  affects  all  attempts  to  measure  and 
interpret  current.  The  interaction  between  the  current  meter's  response  to  the  flow 
and  the  mooring  motion  produce  signals  that  differ  from  the  true  ocean  currents.  For 
lack  of  a  better  term,  we  will  call  the  difference  between  the  signal  and  the  local 
relative  current  "response  error." 

The  response  error  may  arise  from  a  number  of  different  sources,  among  which 

are: 
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1.  nonlinear  response  of  the  flow  sensors  (inertial  effects,  bearing  friction, 
etc.) , 

2.  horizontal  and  vertical  variations  caused  by  instrument  geometry  and  near¬ 
field  flow  patterns  (tie  rod  interference), 

3.  mismatch  of  speed  and  direction  sensors  (vanes/rotors).  For  the  acoustic 
current  meters,  the  primary  response  errors  are  due  to  rotation  of  the  cur¬ 
rent  meter  along  its  axis  due  to  torsional  motion  of  the  mooring  line 
(horizontal  response  error)  and  to  the  tilt  of  the  current  meter  from  a 
vertical  position  (vertical  response  error). 

Three  force  vector  recorders  (FVRs)  were  supplied  by  the  C.S.  Draper  Laborato¬ 
ry,  Cambridge,  Massachusetts  (to  measure  the  tilt  and  torsional  effects  of  the  moor¬ 
ing).  These  FVRs  were  to  have  recorded  acceleration,  magnetic  field,  pressure,  and 
tension  information  at  three  points  on  the  mooring.  The  lower  two  units  failed  to 
yield  useful  data,  and  the  upper  unit  appeared  to  produce  data  of  doubtful  quality 
soon  after  deployment. 

Information  relating  to  the  tilt  of  the  mooring  from  a  vertical  position  was 
deduced  from  the  upper  two  T/P  recorders.  Torsional  motions  could  be  estimated  from 
the  case  direction  measurements  of  the  current  meters.  The  VACMs  returned  a  case  di¬ 
rection  every  sampling  interval  (either  15/16  or  15/8  min),  and  the  acoustic  current 
meters  returned  a  case  direction  every  3  minutes. 

VERTICAL  RESPONSE  ERRORS 

The  vertical  tilt  was  estimated  by  computing 

0  =  arcosine  (aD/L) 

where  aO  is  the  difference  in  the  depth  of  the  upper  two  T/P  recorders  computed  from 
the  observed  pressures,  and  L  is  the  distance  between  the  T/P  recorders  (123  m) .  The 
angles  computed  in  this  manner  are  less  than  10°  for  the  periods  in  which  the  pres¬ 
sure  sensors  were  in  range.  The  error  caused  by  this  10°  tilt  is  about  2% 
(McCullough,  1978),  see  Figure  10;  thus,  the  error  is  expected  to  be  less  than  this. 
Furthermore,  this  error  is  subject  to  very  low  frequency  modulation  and  would  not  be 
expected  to  signif icantly  contaminate  the  high  frequency  region  of  the  spectrum. 

HORIZONTAL  RESPONSE  ERRORS 

Saunders  (1980)  studied  the  horizontal  response  function  of  the  ACM.  Figure  11 
consists  of  plots  of  typical  horizontal  speed  and  response  functions  for  the  ACM.  It 
may  be  seen  that  the  observed  speed  may  vary  between  about  +8%  and  -6%  of  the  true 
speed,  depending  on  the  orientation  of  the  current  meter  case  to  the  mean  flow.  The 
apparent  direction  may  also  vary  between  +5°  of  the  mean  current  direction.  The 
angle  between  the  case  of  an  ACM  and  the  observed  current  is  plotted  as  a  function 
of  time  for  a  typical  case  in  Figure  12.  Long  periods  of  little  case/current  angle 
variation  may  be  seen  interspersed  with  abrupt  large  changes  in  angles.  These 
changes  introduce  the  possibility  of  contamination  of  the  records  by  as  much  as  6-8" 
of  the  mean  speed.  This,  of  course,  is  a  worst  case  situation,  though  one  which 
might  be  expected  to  produce  apparent  high  frequency  contamination  of  the  current 
data  with  spectral  slopes  of  f_2  (step  response).  As  the  existing  data  have  this 
type  of  slope,  it  would  be  difficult  to  detect  this  type  of  contamination  through 
spectral  processing;  however,  it  would  be  expected  to  appear  in  the  records  as  jump 
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CYCLES/HR 


Figure  9.  Rotary  spectra  of  currents  and  mooring  velocities 


Figure  10.  Vertical  response  of  the  Neil  Brown  acoustic  current  meter 
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Figure  1 1 


NBIS  ACM  Horizontal  Response 
Current  Meter  Serial  Number:  08-2296-033 
Offset  angle:  80.50  Degrees 
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Figure  12.  Time  history  of  current  meter  case-current  angle 
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itinuities  in  the  velocity  components.  Some  discontinuities  do  appear,  but  it 
is  not  known  at  this  time  whether  they  are  attributable  to  the  case  response/mooring 
motion  interation. 

To  estimate  the  magnitude  of  the  case  orientation  error,  the  response  function 
[  igure  11  was  applied  to  obtain  estimates  of  the  error  term,  based  on  the  ob- 
erved  currents  at  the  upper  acoustic  current  meter.  The  statistics  for  these  errors 
„  summarized  in  Figure  13,  which  also  depicts  a  "speed-direction"  histogram  of  the 
torsional  error  components.  The  RMS  error  is  1.16  cm/sec  with  a  standard  deviation 
of  0.58  cm/sec.  This  is  about  three  times  as  large  as  the  errors  due  to  translation¬ 
al  motions.  The  magnitude  and  direction  histograms  are  plotted  in  Figure  14.  Figures 
15-17  are  spectral  plots  of  the  observed  currents  at  the  upper  current  meter  com¬ 
pared  with  the  spectra  of  the  angular  response  error.  From  these  plots,  we  can  see 
that  the  torsional  errors  fall  off  with  about  the  same  spectral  slope  as  do  the 
current  components. 

Because  of  the  FVRs  failures,  no  useful  data  that  could  be  used  to  obtain  high 
frequency  accelerations  of  the  mooring  or  tension  records  to  be  used  in  future  moor¬ 
ing  design  and  modeling  were  returned. 

ENVIRONMENTAL  MEASUREMENTS 

The  background  measurements  collected  during  the  ATOM  1979  experiment  fall  into 
two  broad  categories:  shipboard  and  moored  measurements.  The  data  collected  from  the 
ship  provided  information  on  the  vertical  structure  of  the  oceanic  temperatures,  sa¬ 
linity  and  horizontal  component  of  velocity  fields  during  the  deployment,  and  recov¬ 
ery  phases  of  the  experiment.  The  moored  instruments  provided  a  moderately  long  time 
series  of  temperature  and  horizontal  velocity  over  the  period  between  deployment  and 
recovery  (approximately  26  days  at  a  sampling  rate  of  1/min).  The  ship  survey  data 
will  be  discussed  first,  and  is  followed  by  a  discussion  of  the  moored  time  series 
data. 

SHIPBOARD  SURVEYS 

Ihe  primary  purpose  of  the  data  collected  from  the  ship  during  deployment  and 
recovery  was  to  provide  information  on  the  vertical  structure  of  density  and  cur¬ 
rents  near  the  mooring  site,  and  to  provide  a  relationship  between  temperature  and 
i-ncity  that  could  be  used  in  future  analysis  of  the  current/temperature  data. 

Ihe  primary  instrument  used  for  temperature  and  salinity  determination  was  the 
Neil  Brown  CTD.  The  instrument  was  calibrated  prior  to  the  deployment  cruise,  both 
statically  and  dynamically.  CTD,  serial  number  4017  was  used.  The  average  static  er¬ 
rors  for  temperature,  conductivity,  and  salinity  were  -0.0014°C,  -0.00052  mmho/cm 
and  -0.00014o/oo,  respectively.  The  time  constants,  obtained  from  dynamic  drop 
test.',  were  approximately  243  ms  (temperature)  and  59  ms  (conductivity).  The  large 
ciisatch  between  the  temperature  and  conductivity  time  responses  required  matched 
low- pass  filtering  to  avoid  spiking  in  the  computed  salinity  profiles.  The  drop 
rates  were  typically  about  1  m/sec,  thus,  the  data  were  low-pass  filtered  to  give  a 
1  ..  iesolution  (actually  1  dbar  resolution  was  used)  and  the  edited  and  filtered  CTD 
ditj  were  interpolated  and  recorded  at  1  dbar  intervals.  From  these  standard  files 
were  computed  secondary  quantities:  depth,  salinity,  ay,  c>8»  Brunt-Vaisala  fre¬ 
quency  and  sound  speed. 

The  vertical  structure  of  the  horizontal  currents  were  obtained  from  expendable 
shear  probe  drops.  The  horizontal  components  of  current  are  obtained  by  inference 
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Figure  14.  Histogram  plots  of  velocity  errors  due  to  horizontal  current  meter  response 
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CYCLES/HR 


Figure  15.  Rotary  spectra  of  current  and  angular  response  errors 
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Figure  17.  Spectrum  of  V  components  of  current  and  angular  response  errors 
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,  ‘.vired  EMF  gradient.  This  gradient  arises  from  the  passage  of  the  ocean 
,i,  me  earth's  magnetic  field.  The  principle  of  operation  has  been  docu- 
i-if  jrd  et  al .  (1974) . 

of  the  XCP  is  determined  by  an  empirical  formula  relating  the  depth 
■  .  ^  ! :  1  time.  Twenty-five  concurrent  CTD  casts  and  XCP  drops  were  made  to  test 

r  this  method.  Of  these,  twenty-one  had  acceptable  shear  data  (Table 
.its  are  reported  in  Appendix  A.  It  was  found  that  the  standard  equation 
•  ■  ,u curate  enough  for  use  in  calculations  of  Richardson  number  (Ri).  (The  rea¬ 
son  im  mis  lies  in  the  flatness  of  the  velocity  wave  number  spectrum  down  to 
•o  ‘ ,:t  \!  cpm--if  the  depth  accuracy  of  the  XCP  is  not  considerably  better  than 
.)  i  hm:  tne  velocities  being  uncorrelated  over  scales  larger  than  10  m  will  lead 
)  in'  :  eptably  large  errors  in  the  estimation  of  Ri.) 

mprove  the  depth  resolution  of  the  XCPs,  the  initially  computed  depths  were 
cj.vivmd  in  each  case  where  a  concurrent  CTD  cast  was  available.  The  XCP  data  were 
mo!  • -.m.m,lated  to  1  dbar  intervals. 

culations  of  Richardson  number  both  the  XCP  velocity  data  and  the  CTD 
jia  data  were  low-pass  filtered  with  an  8  dbar  cutoff  wave  length.  The 
.  .  numbers  were  computed  from  these  filtered  data.  Typical  profiles  of 

■  :  i  > a i a  frequency,  du/dz,  dv/dz,  and  log  (1/Ri)  are  shown  in  Figures  18  and 


structure  functions  for  the  profiles  of  the  Richardson  number  and 
■m  ))  were  computed,  and  a  typical  example  is  plotted  in  Figure  20.  The  (un- 
ii.ii  ■■■!■  i  zed)  structure  function  for  a  variable,  (x),  is  defined  as 


S $U;x)  =  <0(x)  -  <)>(x  +  s)]2> 


the  ensemble  average  of  X) 
i elated  to  the  covariance  (when  it  exists)  by: 


S$(s;x)  =  <<f>(x)2>  +  < 4>2 (x+C )  >  -2  <  (x )4> (x+c )> 


i  le  of  the  structure  function  over  the  covariance  function  is  that  it  may 
•i  the  variance  and  covariance  function  do  not.  (For  example,  just  consider 

-  i  o  ;  (x)  =  x) . 

uniformity  in  plotting,  we  have  normalized  S<f,U;x)  by 


4<0(x)  -  <4’(x)>]2> 


jnt  -  .*  !  the  averages  by  the  average  over  all  depths.  Thus,  the  structure  func¬ 

tion  nt  interpreted  as  follows:  for  S<f> ( C )  =  0,  the  process  is  correlated  per- 
f  o  ’ '  ,  ,r.  i a y £, ;  for  S$(C)  =  1/2,  the  process  is  uncorrelated  at  lagC;  and  S<f> ( C )  =  -1 
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means  a  perfect  negative  correlation.  The  structure  function  's  used  in  preference 
to  the  correlation  function,  primarily  due  to  its  greater  inherent  stability  in  such 
nonstati onary  fields  as  the  thermal  and  salinity  fields  in  the  upper  ocean. 

Due  to  low-pass  filtering  with  an  8  m  cutoff,  the  structure  function  with  lags 
less  than  8  m  are  not  well  defined.  For  lags  greater  than  8  m,  it  may  be  seen  that 
the  structure  function  for  both  Ri  and  1 og ( ^ /R i )  are  essentially  0.5,  meaning  that 
Richardson  numbers  with  spacings  larger  than  8  m  are  essentially  uncorrelated. 

Wave-number  spectra  for  N,  U2 ,  Vz ,  Ri  and  log  (1/Ri)  have  also  been  com¬ 
puted.  Segment  averaging  was  used  over  all  records  which  had  a  depth  range  greater 
than  512  m  in  which  the  shear  profile  was  not  contaminated  by  the  presence  of  the 
CTD.  These  spectra  are  plotted  in  Figures  21  through  25.  We  can  see  that  for  wave¬ 
lengths  greater  than  the  low-pass  cutoff,  the  spectra  are  essentially  flat.  This  ob¬ 
servation  is  in  agreement  with  the  results  of  Gargett  et  al.  (1980).  The  implica¬ 
tion  is  that  the  shear  field  is  too  uncorrelated  in  the  vertical  (in  the  environment 
encountered  in  ATOM)  for  either  moored  sensor  measurements  with  separations  of  8  m 
(or  so)  or  XCP  probes  used  in  conjunction  with  CTD  profiles  to  determine  the  back¬ 
ground  Richardson  number  field. 

TIME  SERIES  DATA 

Another  major  facet  of  the  ATOM  ’79  experiment  was  the  requirement  to  collect 
current  and  temperature  time  series  from  a  densely  spaced  moored  vertical  array.  The 
purposes  of  these  measurements  were 

1.  to  determine  the  relation,  if  any,  between  the  high  frequency  variability  and 
inertial  frequency  motions, 

2.  to  determine  the  relation,  if  any,  between  the  meteorological  forcing  and  the 
inertial  frequency  motions, 

3.  to  determine  the  degree  to  which  current  components  in  different  frequency 
bands  and  vertical  separations  are  correlated  (particularly  for  the  higher 
frequency  bands) , 

4.  to  determine  the  effect  of  proximity  to  the  Brunt-Vai sa la  frequency  maximum 
on  the  level  of  high  frequency  variability. 

THE  DATA  SET 

In  order  to  test  the  hypotheses,  a  subset  of  the  total  current  and  temperature 
data  was  selected.  This  data  set  consisted  of  20,000  minutes  (approximately  14  days) 
of  ACM  data  from  each  of  the  13  acoustic  current  meters,  starting  at  0000Z,  20  De¬ 
cember  1979  (Day  35  .0000).  Both  components  of  current  and  one  temperature  channel 
were  used.  The  full  length  of  data  from  the  mooring  was  not  used  because  one  current 
meter  failed  (documented  in  Saunders  and  Green,  1980).  The  VACM  data  were  not  usea 
due  to  the  differing  sampling  interval  and  the  different  horizontal  response.  A  5000 
minute  record  from  the  full  data  set  is  plotted  in  Figures  26-30.  The  eastward  ve¬ 
locity  components  are  plotted  in  Figure  26,  the  northward  velocity  components  in 
Figure  27,  the  speeds  in  Figure  28,  the  temperatures  in  Figure  29,  and  the  tempera¬ 
ture  differences  in  Figure  30. 

RELATION  BETWEEN  HIGH  FREQUENCY  ACTIVITY  AND  INERTIAL  MOTION 

One  of  the  hypotheses  of  the  ATOM  experiment  wa  that  inertial  period  motions 
in  the  upper  ocean  would  produce  a  shear  field  that  wc  modulate  the  high  frequen¬ 
cy  internal  wave  intensity.  The  mechanism  for  stimuli  g  such  activity  would  be 
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TABLE  1 


XCP  &  CTO  Correspondence 


XCP  # 

CTD  # 

Temp  Corn. 

Temp  Com. 

H+ 

good  error 

1  (bad) 

5 

X 

2 

- 

X 

3  (bad) 

- 

X 

4 

6 

/ 

/ 

5  (bad) 

7 

X 

6 

7  (up) 

X 

7 

• 

X 

8 

8 

/ 

/ 

9 

9 

/ 

/ 

10  (bad) 

- 

X 

11 

10 

/ 

/ 

12 

11 

/ 

/  ?  anomaly 

13 

12 

>7 

/ 

14 

13 

/ 

/ 

15 

14 

/ 

/ 

16 

15 

/ 

/ 

17 

16 

/ 

/ 

18 

17 

/ 

✓ 

19  (bad) 

18 

X 

20 

18 

/ 

/ 

21 

19 

/ 

/ 

22 

20 

/ 

/ 

23 

21 

✓ 

/ 

24  (bad) 

25 

22 

X 

26 

23 

*7 

/ 

27  (?) 

24 

X 

28  (?) 

25 

X 

29 

26 

/ 

30  (bad) 

27 

X 

31 

27 

✓ 

/ 

32 

28 

/ 

/ 

33 

-- 

X 

/ 

34 

/ 

/ 

35 

31 

/ 

✓ 

36 

32 

/ 

/ 

37 

33 

/ 

/ 

38 

34 

/ 

/ 

39 

35 
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-- 

X 
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Figure  19.  Typical  profile  of  -Log^g  (Richardson  number) 
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Figure  20.  Typical  Richardson  number  structure  functions 
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Figure  21.  Vertical  wavenumber  spectrum  of  Brunt-Vaisala  frequency 
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Figure  22.  Vertical  wavenumber  spectrum  of  U  shear 
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Figure  23.  Vertical  wavenumber  spectrum  of  V  shear 
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Figure  24.  Vertical  wavenumber  spectrum  of  Richardson  number 
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Figure  26.  Eastward  current  components  for  all  acoustic  current  meters  -  first  5000 
minutes 
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Figure  29.  Temperatures  for  all  acoustic  current  meters  -  first  5000  minutes 
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PLOTS  OF  ATI 


2 

4- 

5 

5- 

6 

6- 

7 

7- 

8 

8- 

9 

9- 

10 

10- 

11 

11- 

12 

12- 

13 


o 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 


iW 


7~ ~ 

fv.  WL;  Vam:  '  y. 


A*i 


s*/  V/vvy  1 


\ 


v'A/V^.v/i/i  i '  \i 


•A/, 


”  ’  ,4\  a  /\y  1  ' 


^'^Vv.aV^'v*‘''Av'-\Y' 

^ >yW\ i  V  i*v-' 


-4^- — 

r  *  i  u  *  v  . 


V»( 


rv~, 


^  i  »*y 

jW\>^ ^  v 

Jwy/y/V%|^ 

M,* /V^^Va,/wV^^ 


e._ 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 


0 


TIME 


5000 


Figure  30.  Temperature  differences  between  adjacent  pairs  of  acoustic  current  meters 
first  5000  minutes 
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local  regions  of  intense  shear  that  could  lead  to  dynamic  instability,  overturning 
and  the  eventual  radiation  of  high  frequency  internal  waves  from  the  overturn  re¬ 
gions.  A  second  mechanism  which  could  modulate  high  frequency  internal  waves  is  the 
elastic  scattering  process  (McComas,  1977;  McComas  and  Bretherton,  1977;  McComas  and 
Muller,  1981).  This  is  the  process  by  which  downward  propagating  high  frequency 
waves  are  scattered  by  low  frequency  inertial  waves  with  twice  the  vertical  wave¬ 
length  of  the  high  frequency  waves.  The  induced  diffusion  mechanisms  denote  the 
scattering  of  a  high  vertical  wavenumber,  high  frequency  wave  into  a  nearby  high 
frequency,  :mgh  wavenumber  wave  by  a  low  frequency,  low  wavenumber  wave  (McComas, 
1975;  McComas  and  Muller,  1981).  The  parametric  subharmonic  instability  mechanism 
(McComas,  1975)  generates  low  frequency,  high  vertical  wave  number  waves  from  a  low 
vertical  wave  number  wave  of  higher  frequency.  The  McComas  and  Muller  (1981)  theory 
is  based  on  assumptions  that  the  transfers  of  energy  (action)  are  slow  (occur  over 
many  wave  periods)  and  are  not  influenced  by  steady  vertical  shear.  In  the  upper 
ocean,  particularly  in  the  ATOM  experiment,  the  ambient  steady  shear  is  usually  in¬ 
tense;  consequently,  it  is  doubtful  that  the  weak  interaction  theory  is  directly 
applicable.  Another  assumption  of  the  McComas  and  Muller  theory  is  that  the  spectral 
exchange  processes  are  in  steady  equilibrium.  The  surface  layer  and  the  near  surface 
pycnocline  are  assumed  to  be  source  areas  for  internal  waves;  thus,  generation  of  a 
concentration  of  wave  action  (energy/frequency)  by  a  near  surface  process  momentari¬ 
ly  may  disturb  the  spectral  equilibrium  and  cause  an  exchange  from  high  frequency 
waves  with  high  and  low  wave  number  components  to  lower  frequency.  If  a  process  such 
as  this  occurs,  we  should  expect  to  see  a  strong  correlation  between  high-  and  low- 
frequency  power  levels,  or  stated  as  an  hypothesis:  the  power  levels  of  the  highl¬ 
and  low-frequency,  near  inertial,  internal  gravity  wave  are  correlated. 

The  hypothesis  was  tested  by  correlating  the  energy  envelopes  of  high-passed 
current  data  with  the  envelopes  of  the  energy  of  the  same  currents  band-passed  near 
the  inertial  frequency.  Figure  31  depicts  the  high  frequency  energy  (~1  cph)  envel¬ 
ope  in  the  inertial  band  and  Figure  32  depicts  the  corresponding  inertial  envelopes. 

Referring  to  Figure  31,  it  is  clear  that  the  level  of  energy  in  the  high  fre¬ 
quency  bands  changes  both  with  time  and  depth.  In  general,  the  energy  level  changes 
do  not  appear  to  be  correlated  even  between  pairs  of  adjacent  current  meters.  If  the 
fluctuations  of  the  high  frequency  energy  are  compared  to  those  of  the  inertial  band 
energy,  there  does  not  appear  to  be  any  strong  correlation.  This  is  borne  out  by  ac¬ 
tual  correlation  cal culations ,  the  results  of  which  appear  in  Table  2.  Thus,  we  can 
state  with  some  confidence:  in  the  context  of  the  ATOM  '79  oceanic  environment,  the 
high  frequency  activity  does  not  appear  to  be  strongly  correlated  with  low  frequency 
inertial  wave  activity. 

METEOROLOGICAL  FORCING  OF  INERTIAL  WAVE  ACTIVITY 

The  forcing  of  inertial  wave  activity  by  meteorological  events  (winds)  is  well 
known  (Pollard,  1970;  Pollard  and  Millard,  1970;  Fu,  1981).  The  extent  of  such  forc¬ 
ing  below  the  mixed  layer  in  a  stong  frontal  region  is  not  well  known.  It  was  a  hy¬ 
pothesis  of  this  study  that  the  wind  forcing  at  the  sea  surface  might  be  the  primary 
source  of  inertial  wave  activity  variation  in  the  pycnocline  and  below. 

To  test  this  hypothesis,  the  mooring  was  located  near  (~5  nm)  a  NOAA  Data  Buoy 
Office  data  buoy.  The  winds  recorded  at  this  buoy  and  a  neighboring  buoy  to  the  west 
are  plotted  in  Saunders,  Green,  and  Bergin  (1980).  The  buoy  failed  for  a  period 
during  January,  leaving  a  large  data  gap,  thus  the  wind  forcing  must  be  deduced  from 
the  secondary  data  buoy.  Where  the  data  existed,  they  were  filtered  at  the  inertial 
band,  the  energy  density  computed,  and  the  envelope  computed.  The  energy  envelopes 
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PLOTS  OF  ENERGY 


Figure  32.  Envelope  of  variance  of  inertial  frequency  currents  -  350  hr 


TABLE  2 


1  ■  ■  ■  -  - 

Current 

Meter 

No. 

_ 

HF/INERT* 

INERT/WIND}** 

HF/WINDj** 

1 

.0034 

.534 

.0786 

2 

.662 

.0351 

3 

.603 

4 

.0477 

.608 

.0469 

5 

.0306 

.590 

6 

.0207 

.436 

■  ■ 

7 

.0313 

.497 

.0339 

8 

.0090 

.417 

9 

.0517 

.326 

.0326 

.496 

.0288 

11 

.0458 

.526 

12 

.0373 

.421 

13 

-.0193 

.510 

.0592 

Correlation  coefficients  between  the  high  frequency  (HF)  current  energy,  the 
inertial  current  (INERT)  energy  and  the  wind  energy. 


*  0  LAG 

**  Maxima  over  +^50  hr  Lags. 


for  both  data  buoys  are  plotted  in  Figure  33.  The  upper  plot  is  unreliable  due  to 
missing  data.  Comparison  of  the  lower  plot  with  the  plots  of  inertial  current  energy 
envelopes  show  no  clearly  apparent  correspondence.  A  large  wind  peak  at  about  100 
hrs  was  followed  slightly  later  by  a  significant  current  maximum,  but  the  second  and 
third  wind  energy  peaks  do  not  appear  to  be  associated  with  clear  events  in  the  cur¬ 
rent  records.  However,  the  moderate  correlations  (Table  2)  indicate  a  significant 
relation  between  the  wind  and  the  inertial  motions  in  general. 


Thus,  the  relation  of  wind  forcing  to  inertial  motion  in  strong  baroclinic  cur¬ 
rents  below  the  mixed  layer  appears  significant  in  general,  though  the  gaps  in  the 
wind  data  at  the  mooring  site  reduces  the  strength  of  this  assertion. 

EFFECT  OF  PROXIMITY  OF  THE  8RUNT-VAISALA  FREQUENCY  MAXIMUM  ON  HIGH  FREQUENCY 
INTERNAL  WAVE  ACTIVITY 


The  current  data  were  high-pass  filtered  with  a  cutoff  at  1  cph.  The  filtered 
data  were  squared  and  averaged  over  1000  minute  intervals  and  over  the  entire  20,000 
minute  record  at  each  current  meter.  The  plot  of  the  mean  high  frequency  variance 
versus  depth  (Fig.  34)  shows  a  net  decrease  of  high  frequency  variance  with  depth. 
However,  there  is  wide  variation  between  the  upper  current  meter  (1)  and  the  next 
lower  meter  (2)  and  between  meters  2  and  3.  Current  meter  7  also  appears  to  have 
high  variance  with  respect  to  its  upper  and  lower  neighbors.  Below  current  meter  10 
there  is  a  net  increase  in  high  frequency  current  variance  with  increasing  depth. 

The  peak  of  the  Brunt-Vai  sal  a  frequency  profile  occurs  just  below  100  m  and 
while  there  are  secondary  peaks,  these  are  small  compared  to  the  main  peak.  Although 
the  current  meter  array  was  forced  by  the  current  field  to  undergo  large  vertical 
excursions:  nevertheless,  the  upper  current  meters  were  always  nearer  the  B-V  peak 
than  the  lower  meters.  Thus,  while  on  the  mean  it  is  strictly  correct  to  assert  that 
the  data  show  a  net  decrease  in  high  frequency  activity  with  increasing  distance 
from  the  B-V  peak,  the  large  variations  in  the  high  frequency  variability  tend  to 
mask  this  effect. 

MOORED  VERTICAL  COHERENCE 


The  complex  coherence  function  of  two  randan  processes,  x(t)  and  y(t),  is  prop¬ 
erly  defined  as: 


y  (f) 


Sxy(f) 

/Sxx(f)  Syy(f) 


where  SXy(f)  is  the  exact  cross  spectrum  between  x  and  y  and  Sxx(f)  and  Syy(f) 
are  the  exact  autospectra  of  x  and  y,  respectively.  These  exact  spectra  cannot  be 
obtained  from  natural  oceanographic  data.  The  best  that  can  be  done  is  to  calculate 
an  estimate  of  the  complex  coherence , y  ( f) ,  defined  as: 


;  (f)  = _ <X(UY*(f)> 

/  <X(f)I*(f)  >  <T(f)T*(f)> 

where  X(f)  and  Y(f)  are  the  Fourier  transforms  of  x(t)  and  y(t). 
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(CM2/S2) 


Figure  34.  Average  energy  content  in  high  frequency  band  versus  depth 
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In  actual  practice  the  Fourier  transforms  are  not  really  those  of  x(y)  and 
y(t),  but  of  w(t)  x(t)  and  w(t)  y(t),  where  the  weighting  function  w(t)  is  known  as 
a  window  and  the  averaging  is  done  over  a  sequence  of  Fourier  transforms  computed 
from  a  number  of  time  segments  of  the  data: 


^(f)  =  FFT{wi(t)x(t)| 

Vf)  ■  r  | 


For  the  coherences  shown  in  this  report,  a  Hanning  window  was  used  (Nuttall,  1981). 
The  estimated  complex  coherence  is  usually  expressed  as: 


y(f)  =  |y(f)|ei*(f) 


where  |y(f)|  is  the  magnitude  of  the  coherence  and  the  phase.  I f |  y(f)|  =  1,  the 
signals  are  completely  coherent;  if  |y(f)|  =  0,  the  signals  are  incoherent.  |y(f) 
is  itself  a  random  variable  and  possesses  a  probability  distribution.  The  theoreti¬ 
cal  probability  density  functions  and  cumulative  distribution  functions  for  y(f)  are 
plotted  in  Figures  35  through  37  for  | y ( f ) j  =  0.9,  0.5  and  0,  respectively.  These 
are  based  on  the  theory  presented  in  Carter,  Knapp,  and  Nuttall  (1973)  for  N=37,  the 
number  of  independent  segments  used  in  this  study  and  computed  by  the  algorithms  of 
Barnard  (1981),  and  Lee  (1981). 

A  sample  set  of  current  and  temperature  coherences  appears  in  Figures  38 
through  51.  From  Figure  37  we  see  that  the  80%  confidence  bands  for  zero  coherence 
are  0.055  and  0.249.  Thus,  a  coherence  estimate  must  be  greater  than  0.249  for  there 
to  be  a  probability  of  10%  or  less  that  the  true  coherence  is  zero.  At  the  minimal 
separations  (i.e.,  7  m)  we  see  that  the  velocity  coherence  is  just  at  or  below  the 
30%  upper  confidence  limit  for  all  frequencies  greater  than  1  eph  and  is  definitely 
below  this  limit  for  the  same  frequency  band  and  greater  current  meter  separation. 

The  coherence  functions  between  the  temperatures  at  meters  1  and  6  and  the 
pressure  signal  from  the  upper  T/P  recorder  were  computed.  These  are  displayed  in 
Figures  52  and  53.  The  data  in  the  plots  indicate  that  there  is  essentially  zero 
coherence  between  the  temperature  and  pressure  fluctuations  in  the  frequency  band  of 
1-30  cph. 


1.  currents  in  the  high  frequency  band  (1  cph)  are  generally  incoherent  at 
separations  of  7  m  and  above. 

2.  temperatures  possess  greater  coherence  with  separation  than  the  currents, 
possibly  due  to  internal  wave  motions  that  have  leaked  from  the  band  of 
Brunt-Vaisala  frequencies  where  they  exist  as  free  waves. 
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Probability  Density  Functio 


COHERENCE  801  CONFIDENCE  LIMITS 

Number  of  Segments  Averaged  “  73.  2 


True  Coherence  =  0.  9000 
Upper  Confidence  Limit  =  0.9200 
Lower  Confidence  Limit  =  0.8800 
Bias  =  0.0003 


Estimated  Coherence 


Figure  35.  Probability  density  function  and  cumulative  distribution  for  estimated 
coherence  -  73  segments,  true  coherence  =  0.9 
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COHERENCE  90*  CONFIDENCE  LIMITS 


Number  of  Segments  Averaged  =  73.  0 
True  Coherence  =  0.0000 

Upper  Confidence  Limit  =  0.1800 

Lower  Confidence  Limit  =  0.0400 

Bias  =  0.  1 1 70 


Estimated  Coherence 


Figure  37.  Estimated  coherence  -  73  segments,  true  coherence  =  0.0 
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Figure  38.  Coherence  and  phase  between  U  velocity  components:  meters  1  and  2 


METER  NUMBERS  1  AND  3  NUMBER  OF  FREQUENCIES  256 

VARIABLE  V  NUMBER  OF  SEGMENTS  73 

Figure  39.  Coherence  and  phase  between  U  velocity  components:  meters  1  and  3 
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Figure  40.  Coherence  and  phase  between  U  velocity  components:  meters  1  and  4 
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Figure  42.  Coherence  and  phase  between  U  velocity  components:  meters  1  and  6 
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Figure  43.  Coherence  and  phase  between  V  velocity  components:  meters  1  and  2 
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Figure  44.  Coherence  and  phase  between  V  velocity  components:  meters  1  and  3 
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Figure  45.  Coherence  and  phase  between  V  velocity  components:  meters  1  and  4 
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Figure  46.  Coherence  and  phase  between  V  velocity  components:  meters  1  and  5 
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Figure  47.  Coherence  and  phase  between  V  velocity  components:  meters  1  and  6 
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Figure  48.  Coherence  and  phase  between  temperatures:  meters  1  and  2 
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Figure  49.  Coherence  and  phase  between  temperatures:  meters  1  and  4 
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Figure  51.  Coherence  and  phase  between  temperatures:  meters  1  and  6 
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CONCLUSIONS 


The  scientific  method  is  based  on  the  sequence  of  hypothesis  formulation,  test¬ 
ing  and  reformulation,  usually  with  respect  to  a  simplified  model  of  the  process 
under  study.  The  ATOM  experiment  was  intended  as  an  early  step  in  this  loop  within  a 
framework  of  an  oceanic  internal  wave  continuum  in  a  shear  field  as  the  model. 

We  have  formulated  our  hypotheses  and  tested  them.  In  the  process,  we  have  add¬ 
ed  new  hypotheses  and  have  found  that  some  of  our  original  ones  have  not  been  well 
posed.  The  results  of  the  testing  process  have  provided  input  to  the  next  stage  of 
the  formulation-testing  loop  and  are  summarized  below. 

THE  VERTICAL  DEPENDENCE  OF  HIGH  FREQUENCY  HIGH  WAVE  NUMBER  (HFHWN)  VARIABILITY 

We  hypothesized  that:  HFHWN  variability  is  dependent  on  proximity  to  the  peak 
of  the  buoyancy  profile.  In  the  mean  this  hypothesis  is  true.  The  data  are  erratic 
in  the  sense  that  the  variations  about  the  mean  trend  are  very  large. 

MODULATION  OF  HFHWN  VARIABILITY  BY  INERTIAL  CURRENT  SHEAR 

We  hypothesized  that:  HFHWN  variability  is  modulated  by  inertial-band  waves.  We 
found  no  evidence  that  HFHWN  variability  is  correlated  with  intensity  of  inertial- 
band  modulation  below  the  pycnocline. 

INFLUENCE  OF  LOCAL  METEOROLOGY  ON  HRHWN  VARIABILITY 

Hypothesis:  Low  frequency  (inertial)  variability  is  dependent  on  local  meteoro¬ 
logical  conditions.  This  hypothesis  is  true;  the  variation  of  inertial  band  energy 
below  the  pycnocline  correlate  with  variations  in  surface  wind  intensity.  Above  the 
pycnocline  the  variations  of  the  surface  and  internal  wave  energies  are  known  to  be 
well  correlated  with  surface  wind,  but  there  appears  to  be  little  influence  on  the 
HFHWN  intensity  below  the  pycnocline. 

VERTICAL  COHERENCE  OF  HF  VARIABILITY 

Hypothesis:  HF  variability  is  coherent  over  vertical  separations  greater  than 
5  m  but  less  than  40  m.  This  hypothesis  is  true  for  current  fluctuations.  Coherences 
of  current  fluctuations  above  1  cph  are  not  significant  for  meter  separations 
greater  than  35  m. 

Spectral  analyses  of  the  XCP  data  have  shown  that  the  shear  spectrum  is  white 
down  to  scales  of  8  m.  Thus,  any  future  experiment  should  be  desigred  to  sample  cur¬ 
rents  on  spacings  much  smaller  than  8  m  in  order  to  obtain  good  estimates  of  shear. 

CONTAMINATION  OF  HF  MEASUREMENTS  BY  MOORING  MOTION 

Hypothesis:  HF  measurements  (greater  than  1  cph)  from  a  high  tension,  in-line 
mooring  are  significantly  contaminated  by  mooring  motions.  This  hypothesis  can  be 
rejected  in  the  band  from  inertial  to  about  5  cph,  since  the  mooring  motions  are 
more  than  an  order  of  magnitude  smaller  than  the  recorded  currents.  Above  5  cph  the 
estimated  noise  level  of  the  Doppler  acoustic  tracking  system  becomes  large  enough 
(about  one-half)  to  seriously  impact  the  validity  of  estimates  of  mooring  motion 
effects  on  higher  frequency  current  fluctuations.  This  is  a  technical  problem  in  the 
Doppler  system. 
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Figure  52.  Coherence  and  phase  temperature  at  meter  2  and  pressure  at  the  upper 
T  /  P  recorder 
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Figure  53.  Coherence  and  phase  temperature  at  meter  6  and  pressure  at  the  upper 
T  IP  recorder 


REMOVAL  OF  MOORING  MOTION  FROM  MEASUREMENTS 


Hypothesis:  A  significant  amount  of  recorded  mooring  motion  contamination  of  HF 
internal  wave  records  can  be  removed  (using  tracking  records  or  numerical  models). 
The  test  of  this  hypothesis  had  ambivalent  results  because:  mooring  motions  aid  not 
significantly  "contaminate"  measurements  up  to  about  5  cph.  Down  to  the  noise  limit 
of  the  Doppler  acoustic  tracking  system  it  is  possible  (but  not  particularly 
helpful)  to  subtract  the  vector  time  series  of  the  mooring  motions  from  the  measured 
currents.  Above  that  frequency  the  tracking  data  are  not  sufficient.  A  numerical 
model  of  the  mooring  response  could  conceivably  push  the  correction  to  higher  fre¬ 
quencies,  but  the  model  would  not  be  validated. 

In  addition  to  evaluating  these  hypotheses  we  had  to  consider  some  problems 
that  infringed  on  the  analyses  of  the  measurements. 

IMPORTANT  TECHNICAL  PROBLEMS  TO  BE  ADDRESSED  IN  FUTURE  WORK 

XCP 

We  found  that  the  empirical  fall  distance  formula  for  the  XCP  had  significant 
error,  and  the  empirical  measurements  of  fall  distance  are  subject  to  erratic  scat¬ 
ter  that  cannot  be  fully  rationalized.  It  is  clear  that  introduction  of  a  time  mark¬ 
ing  fall  start  would  decrease  offset  depth  error.  Improved  ventilation  of  the  therm¬ 
istor  would  enhance  empirical  fall  rate/distance  estimates.  We  hypothesized  that  the 
XCP  response  was  sensitive  to  the  electric  field  surrounding  the  CTD.  This  hypothe¬ 
sis  was  found  to  be  true.  The  range  of  vertical  influence  appears  to  be  about  +20  m. 
The  effect  is  not  observed  every  time,  and  it  is  easily  identified  in  proximity 
cases. 

ACM/VACM  AZIMUTHAL  RESPONSE 

As  a  result  of  calibrations  made  after  the  original  experimental  design,  we  had 
hypothesized  that  azimuthal  orientation  of  ACMs  and  VACMs  with  respect  to  the  cur¬ 
rent  could  be  a  source  of  significant  velocity  measurement  error.  This  hypothesis 
was  found  to  be  valid;  in  fact,  the  current  measurement  errors  contributed  by  azi¬ 
muthal  case  orientation  uncertainties  are  greater  than  those  attributed  to  mooring 
motions  by  about  a  factor  of  three.  This  problem  can  be  reduced  by  calibrations  of 
the  azimuthal  responses  of  the  meters,  incorporation  of  these  calibrations  in  the 
conversion  to  current  velocities  and  sampling  of  the  case  orientations  as  rapidly  as 
the  currents.  Alternatively,  the  redesign  of  the  cases  of  the  ACMs  and  VACMs  so  that 
they  have  improved  azimuthal  response  would  reduce  this  error. 

MOORING  DYNAMICS  IN  INTENSE  CURRENTS 

ATOM  was  inadvertent ly  deployed  in  an  excepti onal ly  intense  current,  which 
caused  larger  horizontal  and  vertical  displacements  than  had  been  expected.  If  ATOM- 
like  moorings  are  deployed  in  the  future,  the  certainty  of  keeping  instruments  at 
design  depths  could  be  improved  by  decreasing  the  fluid  dynamic  drag  on  the  elements 
of  the  mooring.  The  dynamical  models  were  not  evaluated  for  conditions  encountered; 
work  should  proceed  on  validating  dynamical  models  with  observations  made  under  a 
broad  range  of  conditions--including  strong  currents. 
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APPENDIX  A 


1.  Derivation  of  the  method  for  obtaining  vertical  velocity  from  the  temperature 
and  pressure  records. 

The  basic  equation  for  temperature  at  a  point  in  the  fluid  is 

~r  +  u-vT  =  NV2T 


or 


if  we  assume 


It  +  “E  *  4 +  4  *  N’2T 


N?n  «  |i. 


Then  we  have 


If  we  consider  a  point  on  the  mooring,  the  observed  velocity  is  given  by 
w0BS  =  w  (due  to  mooring  motion)  +  W  (due  to  water) 

■  WI  +  Ww 

Now  Ww  =  Wqbs  +  Wj 
and  Wj  =  -  If 


wqbs  *  “  If/ lz 


So  the  trye  water  velocity  is  given  by 


w 


3T  /3T  -  1  dp 
3t  /  3z  pg  dt 
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APPENDIX  S:  THE  DERIVATION  OF  VERTICAL  EDDY  VISCOSITIES 


1.  Consider  the  horizontal  momentum  equations;  say  for  u. 

It-  +  U*VU  -fv  = - -  Px  +  vV2U 

3t  pQ 

If  we  denote  an  average  by  an  overbar  and  define 


u  =  u  +  u' , 


The  averaged  equation  is  the  same  except  for  an  extra  term  of  the  form 


u' *vu' 


, ,  3u  1 


3x 


+  V 


,  3U 1 
3X 


+  W 


3u  ' 

3  Z 


If  we  assume  that 


u 


3U ' 
3x  ’ 


« 


w 


3U  ' 

3Z 


and  that 


i  3U 1 
3Z 


82  W  U 


( i .  e . 


u'^1  « 

3Z  ^ 


3Z 


U  W 
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then  we  can  define  an  eddy  viscosity  such  that 


_a 

az 


(wV) 


or 


_a 

az 


Eu 


3u  I 

azj 


=  o 


which  will  hold  whenever 


w  V 


au 

az  =  constant 


The  simplest  choice  is  to  assume  V^u  is  defined  when  the  constant  = 


0  or 
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